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We study the non-compact version of the U(l) gauge-Higgs model in three di- 
q ! mensions for m# = 30GeV. We found that, using this formulation, rather modest 

lattices approach quite well the infinite volume behaviour. The phase transition is 
first order, as expected for this Higgs mass. The latent heat (in units of T c 4 r ) is 
r**- ■ compatible with the predictions of the two-loop effective potential; it is an order 

of magnitude less than the corresponding SU(2) value. The transition tempera- 
ture and < (f*if > in units of the critical temperature are also compatible with 
the perturbative results. 
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1 Introduction 



The main reason for the study of the three-dimensional gauge-Higgs system is 
its relation to the full SU(2) x U(l) Standard model at finite temperature. The 
latter has been studied extensively in recent years in connection with the scenario 
of baryon violation at the electroweak scale during the evolution of the Early 
Universe. 

It is well known that perturbation theory is not reliable for the study of 
such models, because of severe infrared divergencies. One promising approach 
has been to reduce the four- dimensional model at finite temperature to an ef- 
fective model in three dimensions. This can be done if the couplings are small 
and the temperature is much larger than any other mass scale in the theory 
[1, 2, 3]. The parameters of the reduced theory are related to the ones of the 
original model through perturbation theory. The reduced theory has some ad- 
vantages over the original one from the computational point of view [5, 6, 7]. It is 
super-renormalizable and yields transparent relations between the (dimensionful) 
continuous parameters and the lattice ones. Moreover, the number of mass scales 
is drastically reduced: (a) the scale T, present in four dimensions is evidently 
absent, (b) one may also integrate out the temporal component A of the gauge 
field, so its mass scale gT also disappears. Thus there are two mass scales less 
and this reduces substantially the computer time needed to get reliable results. 

The model that has already been studied along these lines [5-8] has been based 
on 577(2) with one complex Higgs doublet. (For work on the same model in asym- 
metric four-dimensional lattices, see [15].) The issues studied have been the order 
and the characteristics (critical temperature, latent heat, surface tension, corre- 
lation lengths) of the phase transition, as well as the reliability of perturbation 
theory deep in the broken phase for several values of the Higgs mass. It is in- 
teresting to see how the above findings are affected by two characteristics of the 
study: the compactness of the gauge group and its non-abelian nature. 

The abelian Higgs model has already been studied both in three and four 
dimensions [14]. The compact U(l) model has been studied on the lattice in [10, 
11]. An interesting aspect of the role of the abelian character of the model would 
be to compare its latent heat against the one of the corresponding SU(2) theory. 

We have chosen to concentrate on the non-compact model. To be exact, only 
the kinetic term for the gauge field is written in the non-compact form; for the 
kinetic term of the scalar field we use the compact formulation. This formalism 
has several advantages: 

• It follows closer the continuum theory, so it should be easier to approach 
the continuum limit. 

• Our results show that with relatively small volumes one gets quite close to 
the thermodynamic limit. 

• The spurious U(l) monopoles, present in the compact formulation will not 
be present any more. 
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• The non-compact version is closer to the exact lattice form of the Landau- 
Ginzburg theory of superconductivity. The variable x of our model, defined 
in the text below, corresponds to k, 2 of the above theory. We recall that 
the the phase transition is of first order for type I superconductors (k < 
4g), corresponding to small Higgs mass; the converse holds for type II 
superconductors. 

We have concentrated on the phase transition line. We have chosen to fix the 
Higgs mass to a fixed value (30 GeV), g to 1/3, m w to 80.6 GeV and study the 
characteristics of the phase transition. 

2 Reduction of the four— dimensional theory to 
three dimensions 

The Lagrangian for the U(l) gauge-Higgs model in four dimensions is well known: 

Ud = \f^F^ + \D^\ 2 + m 2 <p*<p + A(^V) 2 (1) 
The action for this model at finite temperature is: 
S[A^r,x),cp(r,x)) = J* 'tZr J d'x^F^F^ + \D^\ 2 + m 2 <p*<p + X(cp*cp) 2 ), (2) 
where (3 = 1/T. 

If the action is expressed in terms of Fourier components, the mass terms are 
of the type: 

[(2vmT) 2 + (k) 2 ]\A,(n,k)\ 2 (3) 
[(27rnT) 2 + (k) 2 ]\^n,k)\ 2 , (4) 

where n = — oo, . . . , oo. 

At high temperatures T and energy scales less than 2nT the non-static modes 
A^{n 7^ 0, k), ip(n ^ 0, k) are thus suppressed by the factor (27rnT) 2 relative to the 
static A^(n = 0, k) and ip(n = 0, k) modes. The method of dimensional reduction 
consists in integrating out the non-static modes in the action and deriving an 
effective action [2, 3]. 

An important remark is that the mass of the adjoint Higgs field is of order 
gT, which is large compared to g 2 T, the typical scale of the theory. Thus one can 
go on one step further and integrate it out using perturbation theory [5, 6, 7]. 

The effective action may then be written in the form: 

S 3 d e ff[Ai(x), ip 3 {x)] = J d^xl^FijFij + \Dnp 3 \ 2 + m 2 ^3V3 3 + \ 3 (Lp* 3 Lp 3 ) 2 } (5) 

The index 3 in (|5|) denotes the 3D character of the theory. The relations 
between the 4D and 3D parameters are (up to 2 loops ): 
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93 = g 2 (p)T, 



(6) 



m!te) = ^ + i(A 3 + ^)T 

ml = \g\lAT 2 - (9) 

We note that f 2m = — 4c?| + 8A 3 #f - 8A^ and c = -0.348725 [5, 6, 7]. 

The couplings g%, A3 of the three-dimensional theory are renormalization group 
invariant because the theory is supernormalisable. The mass parameter m\ con- 
tains a linear and a logarithmic divergence. 

It is convenient to use the new set of parameters (g%, x, y) rather than the set 
{gh^ m D- x, y are defined by [8]: 

x=^ (10) 

at 

It is evident that x is just proportional to the ratio of the squares of the 
scalar and vector masses; on the other hand, y is related to the temperature. The 
parameters x, y can be expressed in terms of the four-dimensional parameters as 
follows [10, 11]: 

x = l^-^l (12) 
2m 2 w 8tt V ; 



1 1 a/3 

y = t-^. + w^\ x + 



9- 



4g 2 3g 2y 8vr ' 

1,8,2, , y/3g^ 1 
+ 16tt2 1 9 + 3 [X+ Stt }} AnVS 
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1 „ 3\/3 L m„ 

. — (log— + c+-)-^ 
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+ y^(- 4 + 8a; - 8 ^)(J g j 2 + c ) ( 13 ) 
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3 The lattice action 



Discretizing the continuum action (|5|) we get: 

s = AE E ^!+^EE^( if M*)-/(*M i[ + J )] 

a: 0<i<j a 0<i 

+ £[(1 " 2fo - 3^)^*(x)^(x) + p R (<p*(x)ip{x)) 2 ], (14) 

x 

where = A{ Aj(x) - AjA^x), Ui(x) = e iA ^ x \ 

Notice that we use the non-compact version for the gauge field as explained in 
the introduction. The naive continuum limit corresponds to the values: f3 g = oo, 

Ph = 3) Pr = 0. 

The lattice parameters and the (three-dimensional) continuum ones are re- 
lated as follows 0: 

' = (15) 
* " f (16) 
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[(-4 + 8x - 8x 2 )(log6/3 g + 0.09) - 1.1 + 4.6x]. (17) 



16tt 2 

We note that £ = 3.176 at the scale ^ = g\. 



4 The algorithm 

We used the Metropolis algorithm for the updating of both the gauge and the 
Higgs field. It is known that the scalar fields have much longer autocorrelation 
times than the gauge fields. Thus, special care must be taken to increase the 
efficiency of the updating for the Higgs field. We made the following additions to 
the Metropolis updating procedure [8]: 

a) Global radial update: We update the radial part of the Higgs field 
by multiplying it by the same factor at all sites: R(x) — > e^R(x), where £ G 
[— £, e] is randomly chosen. The quantity e is adjusted such that the acceptance 
rate is kept between 0.6 and 0.7. The probability for the updating is P(£) = 
min{l, exp(2V£ — A5(£))} where AS(£) is the change in action, while the 2V^ 
term comes from the change in the measure. 

b) Higgs field overrelaxation: We write the Higgs potential at x in the 
form: 

V(<p(x)) = -a ■ F + R 2 (x) + (3 h (R 2 (x) - l) 2 (18) 

where 
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a _ ( R{x) cosx(x) \ 
y R(x) smx(x) J ' 

F _ ( PhJ2i R{x + i) cos(x(f + + 0(f)) \ 
V PhZiR(x + t)sm( X (x + i) + 6(x)) )' 

We can perform the change of variables: (a, F) — > (X, F, Y) ,where 

F = |F|, f = ; F =, X = a.f, Y = a-Xf. (19) 

The potential may be rewritten in terms of the new variables: 
V(X, F, Y) = -XF + (l + 2^(Y 2 -l))X 2 + Y 2 (l-2/3 fi )+^(X 4 + Y 4 ). (20) 
The updating of Y is done simply by the reflection: 

Y ->• Y' = -Y. (21) 
The updating of X is performed by solving the equation: 

)- expM F', YO) = - exp,-P (X, F, Y)). 

(22) 

The change X — > X' is accepted with probability: -P(X') = min{Pcb 1} ; where 

p _ dV(X,F,Y) / dV(X',F',Y') 

r ° — ax / ax' 



5 Results 

For our Monte-Carlo simulations we used cubic lattices with volumes V = 12 3 , 
16 3 , 24 3 . For each volume we performed 60000 to 110000 thermalization sweeps 
and 70000 to 120000 measurements. We have set the value of x equal to 0.0463. 
According to the relation (|T2|), using mw = 80.6GeV and g — ^, this value 
of x corresponds to a Higgs field mass tuh = 30GeV. We used two values for 
P g , namely (3 g = 4 and /3 fl = 8. For each value of f3h we use the relation QTB|) 
to determine the corresponding (3r. This value of x has been used in references 
[10, 11] in the study of the compact U(l) model, so we use the same value to 
facilitate comparison. The two models should be close for large values of f3 g , 
where the compact formulation probably approaches the non-compact one. The 
phase transition is expected to be of first order, since the mass of the scalar field 
is safely low. 

We used four quantities to locate the phase transition points: 
1. The distribution N(E Hnk ) of E Hnk . 
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2. The susceptibility of Eu n k = J2x,i £l*(x)Ui(x)Q(x+i) (we have set (p(x) 
R{x)jx{*) = R( x )Q(x) ): 

S(Eu nk ) = \/(< (Eu nk ) 2 > — < En n k > 2 ). 



3. The susceptibility of i?2 = ±E X R 2 (%) ■ 

S(R2) = V(< (R2) 2 > - <R2 > 2 ). 



4. The Binder cumulant of Ei in ^: 

n(T? s_ 1 < {ElinkY > 
^K^link) — 1 - q / T-i x 2 ^2' 

3 < (E link y > 2 



We have searched for the (pseudocritical) (3h values yielding (a) equal heights 
of the two peaks of the distribution N(Eu n },), (b) the maxima of the quantities 
S(Eu n },), S(R2) and (c) the minima of the cumulant C(Ei in k). Of course, the 
values PI (A, V) found using each of the above four quantities, depend on the spe- 
cific quantity (denoted by A) which has been employed, as well as on the volume 
V. While searching, we have made use of the Ferrenberg-Swendsen reweighting 
technique [4] to find the pseudocritical f3h for the volume 24 3 . 

In figure 1 we show an example of the distribution of Eu n k in a 16 3 lattice for 
f3 g = 8 and three values of f3h '■ the pseudocritical one (0.337000), one somewhat 
smaller and one somewhat larger than this value. This is just to illustrate the 
way in which the equal height criterion for the critical point has worked. It is 
clear from the figure that the pseudocritical fih yields two maxima of equal height 
in the distribution. For the "small" fth the peak in the region of small values of 
Ei in k is more pronounced, while for the "large" f3h it is the other way around. 
The picture of the two well separated peaks at criticality is the signature of a first 
order phase transition; the two peaks correspond to the two coexisting metastable 
states. 

In figures 2 and 3 we depict the behaviour of the susceptibilities S(Eu n k) 
and S(R2) for j3 g = 4 versus j3h for three lattice volumes. We have fitted curves 
through the data and show them in the figures; for the 16 3 and 24 3 lattice volumes 
we also give the actual measurements. It is evident that the curves represent the 
data quite nicely. To calculate the error bars we first found the integrated auto- 
correlation times Ti n t(A) for the relevant quantities A and constructed samples of 
data separated by a number of steps greater than Ti nt (A). The errors have been 
calculated by the Jackknife method, using the samples constructed according to 
the procedure just described. We observe that the peak values for the suscepti- 
bilities increase almost linearly with the volume in both cases, which is evidence 
for a first order phase transition. In figure 4 we depict the behaviour of the 
Binder cumulant C(Eu n k) at (3 g = 8 for three lattice sizes. We again show the 
real measurements for 16 3 and 24 3 only and just give the fitted curves for 12 3 . The 
error bars have been calculated by the Jackknife method [13], in the same way 
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Figure 1: Distribution of Eu n k 
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Figure 2: Susceptibility of Eu n k 
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Figure 3: Susceptibility of R2 
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Figure 4: Binder cumulant of Eu n k 
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as in the case of the susceptibilities. The volume dependence of the cumulants is 
again characteristic of a first order phase transition. 

The use of finite lattices is the reason why the /3^(A, V) values that we have 
found employing the various criteria are slightly different. Thus, one should 
extrapolate these values to infinite volume. We have adopted the ansatz: 

P*h(A V) = W(oo) + ^, 

The constant c(A) is expected to depend on the quantity A, while the extrapo- 
lated value P^ r (oo) should not depend on A; that is, the infinite volume extrap- 
olation for the critical point should not depend on the quantity used. 

Figures 5 and 6 deal with the extrapolation to infinite volume for j3 g = 4 and 
P g = 8 respectively. They contain the data for the pseudocritical /3£(A, V) values 
obtained from the various quantities A versus the inverse lattice volume, along 
with the linear fits to the data. The error bars in (3h have been found from the 
statistical error of the values of the quantities A at the critical point. We note 
in passing that, at finite volumes, the smallest pseudocritical values are given by 
the cumulant of Eu n k, then follow, in ascending order, the values given from the 
equal height, the susceptibility of Ei in k and the susceptibility of i?2.This holds for 
both values of f3 g . One may also observe that the infinite volume extrapolation is 
almost independent from the specific quantity used: the differences at the point 
y = between the various extrapolated critical values are less than 1CT 5 . To 
be specific, the critical values lie in the interval (0.340295, 0.340298) for f3 g = 4 
and in (0.336932, 0.336940) for (3 g = 8. In reference [10], treating the compact 
U(l) model, the best pseudocritical value correesponding to our results has been 
obtained for (3 g = 8 for a 32 3 lattice. From the relevant figure one may read out 
a pseudocritical value about 0.3370, a result consistent with ours. These results 
suggest that the non-compact formulation allows one to obtain similar results to 
the ones of the compact formulation in a quite economical way. 

The next quantity we are going to deal with is the critical temperature. It 
may be determined by noting that, for each /3 g , the quantity /3£ r (oo) yields y cr 
through equations Q16] , |T7D ; then equation ([13]) gives T cr . The results are to be 
found in table 1. We mention that in reference [8] the critical temperature for 
the SU(2) model at (3 g = 8 and m H = 35GeV has been found 9A.181GeV. No 
lattice result is reported for this quantity in the paper ]ll| on compact U(l), but 
there is the result 148.83GeV at m# = 35GeV from the perturbative effective 
potential. We will say more about this later on, but we remark at this point that 
this value is quite close to ours. 

Having found the critical temperature we estimated the latent heat, that is 
the energy released in the transition. We have used the formula [8, 9]: 

T 1 /If 2 

7T = ~-7ft&^A < *2 > . (23) 

cr cr 

We note that g\ = g 2 T cr . The quantity A < R2 > is the difference of the 
R2 expectation values between the phases. We have measured the values of 



10 



0.34045 



0.3404 



0.34035 



0.3403 



0.34025 



Mhiggs=30GeV Bg=4.0 



0.3402 



susceptibility R2 hs — i 

susceptibility Elink h — i 

equal height Elink i-b— i 

cumulant Elink hx — i 



0.0001 0.0002 0.0003 0.0004 0.0005 0.0006 0.0007 0.0008 

1/V 



Figure 5: Extrapolation for (3 g = 4 
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Figure 6: Extrapolation for (3 g = 8 
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Figure 7: Extrapolation for A < R2 > 
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Figure 8: Extrapolation for A < R2 > 



A < R2 > from the R2 distributions for each lattice volume at the three different 
pseudocritical values of /3^(A, V) [12] the quantities A being the susceptibilities of 
Ei in k and R2 and the equal height signal of N(Eu n k). Figures 7 and 8 show these 
sets of three measurements versus the inverse volume of the lattice for the cases 
P g = 4 and P g = 8. The error bars are due to the uncertainty of each pseudocritical 
value, as well as to statistical dispersion; they turn out to be rather big, especially 
for the smallest volume. The final values of A < R2 >, which have been used for 
the calculation of have been found from a linear fit of these data with ^. The 
convergence of the three straight lines to a common value in the limit of infinite 
volume is fairly good. The results for the quantity ^ can be found in table 1. 

-*- cr 

The corresponding results for the SU(2) case at ma = 35GeV have been reported 
to be 0.180 ± 0.001 in 



121 and 0.256 ± 0.008 in 



The former result has been 
obtained for P g = 12, while for the latter an extrapolation to P g = oo has been 
done. We observe that our values of this quantity for U(l) are less than the one 
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Figure 9: 2-loop effective potential. We use phi = (pyT 



tenth of the values for SU(2). This permits one to be sure that the U(l) part of 
the Standard Model gauge group plays only a secondary role in the scenario of 
the Electroweak Phase Transition. The relative factor of ten is so big, that this 
conclusion cannot be spoiled by the rather large errors. 



Table 1 


P B 


T 




<(ff(p> /T cr 


4 


131.50(3) 


0.0135(4) 


0.255(6) 


8 


131.18(14) 


0.0172(6) 


0.308(10) 


pert. 


132.64 


0.0150 


0.285 



For the (3-dimensional) lattice quantity < {p*{p > /T cr , also appearing in table 
1, it is important that one subtracts the "infinities" from the lattice results. The 
relevant formula reads 101: 



<(p*(p>= -(3 h (3 9 gt <R2> 



4tc 



8tt 2 



log6(3 g + C + y~ 6 



where £ = 3.176, C = 0.09 and 5 = 1.94. 

One may use the U(l) effective potential || to determine the critical temper- 
ature T cr , as well as the quantity < ip*ip > /T cr and T^t and compare with the 

•* cr 

corresponding quantities from the lattice. (We note that the critical tempera- 
ture is defined in perturbation theory by the equality of the two minima of the 
potential.) The perturbative predictions are also displayed in table 1. 
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Figure 9 depicts the two-loop effective potential versus the four-dimensional 
scalar field (a) for the critical temperature T cr>pert , and (b) for two other neigh- 
bouring temperatures, one corresponding to the symmetric phase and the other 
to the broken phase. 

In principle one should perform the extrapolation to large values of [3 g . How- 
ever one is not sure about the exact (3 g dependence of the various quantities, so 
we postpone this until we get results for even bigger j3 g . 

We observe that the lattice T cr is smaller than the prediction from pertur- 
bation theory and decreases with f3 g (in agreement with the SU(2) results for 
small mu [8]). The other two quantities that we measured, namely ^4- and 

-*- cr 

< tp*(p > /To-, increase with f3 g and their values are compatible with the pertur- 
bative ones. 
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